Vortex motion in Nb/PdNi/Nb trilayers: new aspects in the flux flow state 
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Abstract 

We study the dynamics of vortex lines in Supercondutor/Ferromagnet/Superconductor (SFS) heterostructures at mi- 
crowave frequencies. We have employed swept-frequency, Corbino-disk and resonant, dielectric-resonator techniques 
to obtain the field and temperature dependence of the vortex-state parameters. We concentrate here on the genuine 
flux-flow resistivity p/y , that we access at subcritical currents using a sufliciently high driving frequency. We find that 
Pff does not follow the well-known Bardeen-Stephen model. Instead, it is well described by a full time-dependent 
Ginzburg-Landau expression at very thin F layer thickness, but changes to a previously unreported field-dependence 
when the F layer exceeds a few nm. 
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1. Introduction 

The dynamic properties of the vortex state have been 
studied extensively with experiments, theory and nu- 
merical simulations, with renewed interest after the dis- 
covery of cuprate superconductors. Complex vortex 
matters have been identified, and their properties the- 
oretically investigated [ilj. In face of the enormous 
amount of progress in the field, the simplest, free flux 
flow (FEE) has received much less attention, despite its 
direct connection to the physics in the vortex core. The 
discovery of Iron-based superconductors, and the re- 
vamped interest in superconductor/ferromagnet hybrid 
structures make this issue even more interesting than be- 
fore. 

Experiments on the vortex motion resistivity often 
compare the flux-flow state to the prediction of the 
Bardeen-Stephen (BS) model, which predicts 

Pff,BS = cp„B/iuoHc2 (1) 

where p„ is the normal state resistivity, B is the field in- 
duction, Hc2 is the temperature-dependent upper critical 
field and c ~ 1. This behavior is often observed, apart 
the field region approaching Hc2, where additional ef- 
fects, e.g. the suppression of the condensate, increase 
Pff Sihowe. p„B/i2oHc2- 
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In this paper we study the flux flow resistivity pff 
in SFS trilayers using high-frequencies currents, of the 
order of several GHz. The motivation resides in the pos- 
sibility to access the genuine flux-flow state (defined as 
irrelevance of pinning and free motion of vortices, hin- 
dered solely by the viscous force) without the need for 
large currents. In fact, motion of vortices is a complex 
matter that involves interaction with defects, relation 
with the vortex lattice elastic moduli, and the nature of 
the vortex core (responsible for the viscous drag). How- 
ever, high frequency measurements prove useful in sim- 
plifying the physical problem: at several GHz vortices 
oscillate over distances that are typically much shorter 
than the vortex lattice constant and of the defect average 
distance |2 |. Thus, they remain mostly near their equi- 
librium position. In most cases one can apply safely 
the mean-field theory developed for vortex dynamics 
El m 13] El . Interestingly, a single formulation is able 
to describe many of the mentioned theories. In fact, the 
vortex motion resistivity can be described by the follow- 
ing expression IT): 

E + Hv/v) 

Pvm - Pvm,\ + lPrai,2 - PffZ . , (.^J 

\ + i(v/v) 

where the resistivity pff represents the free flux flow 
value, reached in the high frequency limit in which vor- 
tices experiences only the dissipative viscous drag, v 
is a characteristic frequency and the dimensionless pa- 
rameter < e < 1 is a measure of thermal activation 
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phenomena. For the purposes of this paper, we expUc- 
itly note that, when the measuring frequency v » v, 
Pvm - Pff, irrespective of pinning and creep (whose 
relevance are within v and e). Since in most supercon- 
ductors the characteristic frequencies lay at or below the 
microwave frequencies, an experimental investigation at 
several GHz is a powerful tool to access the genuine p//. 

2. Samples and experimental setup 

Nb/Pd84Nii6/Nb trilayers have been grown by UHV dc 
magnetron sputtering at room temperature, as reported 
elsewhere |[8|. Samples were grown onto sapphire sub- 
strates to avoid the deleterious contribution of semi- 
conducting substrates to the microwave measurements 
||9l . The nominal thickness, calibrated with the depo- 
sition rate, was as follows: all trilayers had Nb thick- 
ness d^b-^^ nm for both upper and lower layers. In 
this work we concentrate on the samples with PdNi 
thickness dp- 1, 2, 9 nm. dc resistivity measurements 
in samples grown on the same conditions consistently 
yielded a constant normal state p„ ^22 + 5 yuQcm for 
temperatures < T < 30 K. Critical temperatures of 
6.25 K, 5.1 K and 4.25 K for the samples with dp- 1, 
2, 9 nm, respectively, were evaluated by the disappear- 
ance of the microwave signal and were consistent with 
dc measurements with a criterion of 90% of the normal 
state. 

We have used two different setups in order to mea- 
sure the microwave response of SFS heterostructures. 
The Corbino disk setup was used to obtain p fj from the 
frequency dependence of p„„ between 2 and 20 GHz, 
in SFS structures with dp -I and 2 nm. The setup 
and the measuring technique has been extensively de- 
scribed previously [lO], and we report here only the 
results. In this context, it is important to stress that 
the swept-frequency data directly yielded a character- 
istic frequency v which decreased with dp, and was as 
low as 4 GHz in the sample with dp -2 nm. Thus, we 
have chosen to perform measurements in trilayers with 
dp - 2 and 9 nm with a more sensitive dielectric res- 
onator at the operating frequency vq - 8 GHz. 

The dielectric resonator (DR) is based upon a rutile 
(Ti02) cylinder, of approximate dimensions 2.73 mm 
(height) X 3.9 mm (dia.), enclosed in a OFHC Cu cylin- 
drical enclosure. The dielectric rod is placed onto the 
sample, which in turn is partially masked by a thin metal 
mask with a circular hole in order to maintain the az- 
imuthal symmetry, and it is in close contact with the 
bottom wall of the enclosure. The high permittivity of 
rutile {er ~100, for a detailed characterization see UJJ) 
ensures that the evanescent electromagnetic field drops 
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Figure 1 : (a) sketch of the section of the DR resonator (proportions 
are approximate). Rutile rod: hatched area; sample: black rectangle, 
(b) log plot of the resonance curve for the trilayer with dp =2 nm at 
T =4.08 K and in zero field, (c) field-variation of the resonance curve. 



quickly to zero out of the dielectric. The DR works in 
transmission with very weak coupling, thus making cor- 
rections of the directly measured quality factor Q unnec- 
essary fl2|. The measurements here reported are all in 
the linear regime, as directly checked by us by varying 
the incident power over a factor of 10 without observing 
any variation in the response. A sketch of the resonator 
is reported in Fig{T^. The resonator was placed in a He- 
flow superconducting cryomagnet, and magnetic fields 
up and over the upper critical field Hc2 could be applied 
perpendicular to the sample plane. At fixed tempera- 
tures (stability < 0.001 K) the field was slowly swept. 
At each magnetic field the resonance curve was mea- 
sured with a Vector Network Analyzer and fitted to a 
Lorentzian, yielding Q and the resonance frequency vq. 

All samples were thin enough to ensure the realiza- 
tion of the thin-film approximation in the microwave 
measurements 1 13 1. Thus, the field variation of the com- 
plex resistivity Ap{H) = piH) - p(Q) = Api -i- iAp2 was 
obtained from Q and vq as: 
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where G is a geometrical factor and d = 2dNb + dp is the 
total thickness of the sample. The reduced thickness im- 
plies that screening is vanishing small approaching T^, 
so that data close to the normal state cannot be measured 
with the present technique. We note that G and d act as 
mere scaling factors. In Figure [T] we report the evolu- 
tion of the resonance curve upon sweeping the magnetic 
field. Q and vq change with the field (the resonance 
curve broadens and the peak shifts), reflecting changes 
in pi and p2- From the application of Eq.s Q, ([3]l and 
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the procedure described in |7 1 we found that, to a good 
approximation, pff ^ Ap\ (being Api a lower limit for 
Pff when the correction could be detected). For the pur- 
poses of this paper, we could safely take p// - Ap;. 

In some samples and in a Nb sample (for reference) 
we checked at several temperatures that the microwave 
resistivity as measured with the DR and the Corbino 
disk yielded the same results. 
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Figure 2: Normalized flux flow resistivity vs. reducedfield, py^//p„ vs. 
H/Hc2, at approximately the same reduced temperature T/T^ k 0.8 
for tliree trilayers witli df = 1,2, 9nm. Continuous line: BS predic- 
tion, Eq.l[TJ, with c = 1. None of the samples follows the BS model. 
Inset: enlargement at small reduced fields. 



3. Results and discussion 

In Fig|2]we report the normalized flux flow resistivity 
in the three samples studied at similar reduced temper- 
ature T/Tc ~ 0.8, as a function of the reduced magnetic 
field. The sample with df =1 nm has been measured 
with the Corbino disk. In this plot in reduced units one 
would expect a similar, if not the same, behavior for 
Pff. However, despite the same reduced temperature, 
the three samples behave dramatically difi'erent. It is 
apparent that, with increasing df, the field dependence 
of Pff progressively changes curvature, from upward 
(df -I nm) to downward (df = 9nm). Even more strik- 
ing, the large diff'erence between the samples with df - 
1 and 2 nm: a close analysis (see the enlargement in the 
inset of Fig. |2]) shows that even at low fields there is no 
proportionality py^y oc H, thus invalidating the BS model 
in all the samples. 

We discuss first the sample at df =1 nm. First of 
all, it should be made clear that the upward curvature, 
and the fact that pff < H/Hc2 cannot be ascribed to 
pinning or creep: the present measurements have been 
taken with the Corbino disk, and the full frequency de- 
pendence of the real part of the resistivity has been 



recorded (typical frequency sweeps, and fits by Eq.([2]i, 
have been reported in [10|). The upward curvature re- 
flects an overestimate of pff by the BS model. We 
then compared the results to more refined theories for 
the FF. We compared the data with the microscopic the- 
ory by Larkin and Ovchinnikov (LO) 1 15 1, but were un- 
able to obtain more than qualitative agreement in the 
low field and H — > Hc2 limit. However, we found that 
a full calculation within the time-dependent Ginzburg- 
Landau framework (TDGL) lfT6l matched very well our 
data. The calculation yielded, for the normalized flux 
flow resistivity: 

PJl . 1 (4) 

p„ 1 + (i^oHci - B) laB 

where a » 0.4, as calculated in IfTTl . The BS model is 
recovered for a - 1 . In Fig j3]we report the same data for 
the sample with df - \ nm as in Figp] together with the 
curves calculated on the basis of Eq^Q for a =0.3, 0.4, 
0.5. To check further the agreement between the data 
and the TDGL expression, we measured pyy at a differ- 
ent temperature. As can be seen, with the simple change 
of the value of Hc2, Eq.(|4| fits the data at different tem- 
perature as well. While we can state that we found 




Figure 3: Normalized flux flow resistivity vs. applied field, Pfflpn 
vs. H, in the sample with df =\ nm, at 5 K (full dots, same data 
as in Fig|2| and 3.58 K (open diamonds). Continuous line: TDGL 
calculation, Eq.jiJ, with a = 0.4. By varying only the value of 7/^2 the 
TDGL calculation fits also data at a different temperature. Short-dash 
and long-dash lines exemplify the sensitivity of the TDGL calculation 
to the parameter a. 

agreement between our pff data and a "conventional" 
theory, so that at least in this respect the trilayer with 
very short df behaves essentially like an ordinary su- 
perconductor, the same is not true for the trilayers with 
larger df. We argue that the increasing df induces a 
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progressive crossover between a "conventional" behav- 
ior and a novel flux flow regime. The data in the sample 
with df -9 nm exceed by far the BS expectation, in- 
dicating that additional losses are induced by the mag- 
netic field. The emerging of this feature can be observed 
at low fields in the sample with dp -2 nm (see inset of 
Fig|2]l. While we have not yet an explanation of the phe- 
nomenon, we can speculate over a possible scenario. At 
very small dp the F layer acts simply to reduce (with 
respect to a Nb sample 30 nm thick, where -7.5 K 
IfTS l) but leaving all the main superconducting features 
unaffected, including the spatial uniformity of the or- 
der parameter. With increasing dp, the order parame- 
ter acquires a modulation across the trilayer (connected 
to the observed oscillations of with dp LI 9. .20, 21 J ). 
The magnetic field may trigger this delicate state, induc- 
ing a relatively wide region where the order parameter 
is substantially depressed. The electromagnetic field, 
uniform along the thickness of thin films, would give 
a response as if part of the trilayer would be essentially 
normal (effective thickness smaller than the geometrical 
thickness), and dissipation larger than expected would 
be observed. While this is at present only a specula- 
tion, it indicates that an extension of the measurements 
in samples with different dp, and possibly with different 
F materials, would be useful to clarify the anomalous 
flux-flow behavior here observed. 

4. Summary 

We have studied the flux flow resistivity of 
Nb/PdNi/Nb trilayers. We have found a conventional 
behavior of pyy, different from the Bardeem-Stephen 
model but well described by a TDGL calculation, when 
the ferromagnetic thickness is extremely small. When 
the thickness of the ferromagnetic layer increases, pjf 
exhibits a previously unreported field dependence. A 
peculiar feature is that it exceeds the Bardeen-Stephen 
value, p„H/Hc2, in almost the entire field range. A 
quantitative explanation of the phenomenon is not given 
at present. 

This work has been partially supported by an Ital- 
ian MIUR-PRIN 2007 project. 
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